library("kernlab")
// data("BostonHousing")



mydata = BostonHousing

# 計算取10%的訓練資料筆數
np = ceiling(0.1*nrow(mydata))
np

# 取出10%為測試資料
test.index = sample(1:nrow(mydata),np)
test.index

my.testdata = mydata[test.index,]	#測試資料
my.traindata = mydata[-test.index,]     #訓練資料


//mymodel <- ksvm(crim ~ ., data = my.traindata, type = "C-bsvc", kernel = "rbfdot", kpar = list(sigma = 0.1), C = 10, prob.model = TRUE)


set.seed(123)
x <- rbind(matrix(rnorm(120), , 2), matrix(rnorm(120, mean = 3), , 2))
y <- matrix(c(rep(1, 60), rep(-1, 60)))
svp <- ksvm(x, y, type = "C-svc")
plot(svp, data = x)


rvmm <- rvm(x, y, kernel = "rbfdot", kpar = list(sigma = 0.1))
rvmm